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It has been widely believed that variability of the fine-structure constant a would imply detectable 
violations of the weak equivalence principle. This belief is not justified in general. It is put to rest 
here in the context of the general framework for a variability [J. D. Bekenstein, Phys. Rev. D 
25, 1527 (1982)] in which the exponent of a scalar field plays the role of the permittivity and 
inverse permeability of the vacuum. The coupling of particles to the scalar field is necessarily such 
that the anomalous force acting on a charged particle by virtue of its mass's dependence on the 
scalar field is cancelled by terms modifying the usual Coulomb force. As a consequence a particle's 
acceleration in external fields depends only on its charge to mass ratio, in accordance with the 
principle. And the center of mass acceleration of a composite object can be proved to be independent 
of the object's internal constitution, as the weak equivalence principle requires. Likewise the widely 
employed assumption that the Coulomb energy of matter is the principal source of the scalar field 
proves wrong; Coulomb energy effectively cancels out in the continuum description of the scalar 
field's dynamics. This cancellation resolves a cosmological conundrum: with Coulomb energy as the 
source of the scalar field, the framework would predict a decrease of a with cosmological expansion, 
whereas an increase is claimed to be observed. Because of the said cancellation, magnetic energy 
of cosmological baryonic matter is the main source of the scalar field. Consequently the expansion 
is accompanied by an increase in a; for reasonable values of the framework's sole parameter, this 
occurs at a rate consistent with the observers' claims. 



I. INTRODUCTION 



Observations of fine structure multiplet splittings in the absorption line systems of distant quasars have lately 
suggested [1] that the fine-structure constant a at cosmological epochs with redshift 1-3.5 was lower than it is today. 
This evidently revolutionary claim begs for a theoretical framework to enable judgment as to whether such a variation 
can be consistent with accepted physical principles. Variation of a was first considered theoretically by Jordan [2], 
Teller [3] and Stanyukovich [4]. 

Already before Gamow's influential speculation that a varies linearly with cosmological time [5], Dicke had made 
the point that any variation of a can be regarded equally well as due to variation of particle charge, or alternatively 
of the speed of light c or of H, with the choice being a matter of convenience, not physics [6]. He also exhibited a 
theory in which the Maxwell invariant couples linearly to a standard massless scalar field as an example of a variable a 
theory. In this theory the permittivity and reciprocal permeability of the vacuum vary in consonance with the scalar 
field; it can also be regarded as a variable charge theory. Independently of the specific theory of a variation, Dicke 
provided an argument [6, 7] that spatial variation of a, which would be expected to accompany cosmological temporal 
variation if the underlying theory is covariant, contradicts the weak equivalence principle (WEP). The essence of the 
argument is that a non-negligible fraction £ of the mass M of any chunk of ordinary neutral matter is Coulombic in 
origin. The Coulomb energy Eq — C^Mc 2 should scale with the square of the constituent charges, meaning it should 
be proportional to a, and should thus depend on position if a does. The force on mass M should thus contain in 
addition to the Newtonian part — AfV^jv an anomalous portion — VE'c = —(dEc/da)Va = — ((Mc 2 /a)Va. Dicke 
conjectured [6, 7] that Vcv/cv w ac~ 2 V '(/)n ■ In this view - here designated the Coulomb model regardless of the 
law of a variation assumed - the acceleration of M comprises an anomalous fraction £cv which varies from material 
to material (for example, the nuclear contribution to £ ranges from 0.0016 for aluminum to 0.0041 for lead), thus 
engendering a violation of the principle of universality of free fall of neutral matter, a special case of the WEP. Uzan [8] 
has given a masterly review of this subject, as well as of the whole question of a variability. 

I made use of the Coulomb model for the anomalous force when drawing conclusions from the general field-theoretic 
framework of a variability I formulated two decades ago [9], The model has also been uncritically adopted by most 
subsequent investigations in the subject [10, 11, 12, 13, 14, 15, 16, 17, 18]; a refreshing exception is Landau et al [19]. 
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Within the general framework the Coulomb model predicts that the anomalous acceleration is a fraction ~ £ 2 of 
the total one; this is within an order of magnitude of the fraction £a implied by Dickc's conjecture (see Ref. [9] 
and Sec. VA below). Similar results are in evidence in other treatments. They have given rise to the widespread 
belief that a variability necessarily implies violations of the WEP detectable by Eotvbs-Dicke-Braginsky (EDB) type 
experiments [8], particularly if the claimed cosmological variability [1] is essentially correct. 

But as shown below, such sweeping conclusion is unwarranted; it all depends on the structure of the underlying field 
theory. Contrary to intuition and my original supposition, in the general framework [9] Coulomb energy of matter is 
found to be unimportant as a source of the scalar field responsible for spatial a variability. Spatial gradients of that 
field are necessarily much smaller than has been generally appreciated heretofore. This finding also impinges upon 
cosmology as it leads to a modification of the widely used cosmological equation for a variability. 

The paper is designed as follows. In Sec. II I recapitulate the general variable a framework and discuss its relation 
to other theories of a variability. I redcrive the equations of motion for scalar and electromagnetic fields and particles 
within the framework, and remark that one is prevented from assuming — as widely done — that the dependence of 
particle masses on the scalar field, which enters both in the anomalous force on particles and in the source of scalar 
field, comes from Coulomb energy. Rather the equations of the theory themselves, being nonlinear, determine the 
nature of their sources. 

This is indeed seen from the exact solution — presented in Sec. Ill — for the fields of a pointlike charge held 
initially at rest in a uniform external electric field. The dependence of the particle's mass on scalar field is not of the 
sort expected from the Coulomb energy model: the anomalous force is cancelled by a correction to the usual electric 
force. Thus there is no indication that different types of particle with like charge-to-mass ratio would move differently 
in the same external fields, and so there is no basis for a violation of the WEP. I also show that within the framework 
a charged particle has a minimum possible extension, and that gravitational corrections to the mentioned results are 
small, even for the most compact charge. 

In Sec. IV I generalize the said solution to the case of many charges held fixed in space. I prove rigorously that 
at an initial moment, the center-of-mass acceleration of a collection of charges (some may be zero) starting at rest 
depends exclusively on its total mass and charge, and on the external electric field, but not on the structure of the 
collection. Again no violation of the WEP is in evidence. 

Sec. V shifts the focus from the microscopic description of individual charges to the macroscopic description in terms 
of smoothed electromagnetic and scalar fields and their smoothed sources. The macroscopic scalar field is shown by 
a variety of examples to be too small to effect WEP violations detectable in the foreseeable future. Ordinary matter 
contributes two terms of Coulomb origin to the source of the scalar field. Taken together uncritically these would 
suggest that in the vicinity of a chunk of matter, Va/a ~ (c~ 2 V4>n, of the same order as Dicke's conjectured spatial 
a variability. But when the mentioned exact solution is taken into account, the two Coulomb terms are seen to cancel 
each other. The source retains only terms of higher order; these are incapable of generating WEP violating signals 
observable in the foreseeable future. 

Casting about for other sources of a variability, I show in Sec. VI that contributions of spin and orbital magnetic 
dipoles in ordinary matter to the sources of the scalar field are likewise too weak to engender violations of the WEP 
at soon-to-be observable levels. 

In Sec. VII I show that the cancellation of the Coulomb energy in the source of the scalar field makes baryonic 
magnetic energy the dominant source of cosmological a variability. This, together with the significant value of the 
fundamental scale of the theory permitted by WEP tests (as newly understood), makes it possible to understand the 
observed cosmological growth of a and the rate of it as reflecting standard properties of cosmological baryonic matter; 
there is no need for this to postulate dark matter with peculiar electromagnetic properties. Sec. VIII summarizes the 
conclusions and caveats on them. 

Below I shall employ the signature {—1,1,1,1} and the convention that Greek indices range from to 3, while 
Latin indices take on values from 1 to 3. The time coordinate is denoted by t or x a . 

II. FRAMEWORK FOR a VARIABILITY 

A. Modifying electrodynamics 

The general field-theoretic framework for a variability [9] is based on eight assumptions: (1) for constant a the 
framework's electromagnetism reduces to Maxwell's with a minimal coupling to charged matter, (2) a dynamics comes 
from an action, (3) this as well as electrodynamics' action are relativistic invariants, (4) the overall action respects 
gauge invariance, (5) electromagnetism is causal and (6) respects time reversal-invariance, (7) any length scale in the 
theory is not smaller than Planck's length lp = (ftG/c 3 ) 1 / 2 « 1.616 x 10 -33 cm, and (8) gravitation is governed by 
the Einstcin-Hilbcrt action. 
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The choice of units that makes G, h and c constant shifts the burden of variation onto the charges. Simplest is 
the case where all charges vary in unison: = e ai e(x A '), where denotes the coupling constant of particle i and 
e(x^) is a dimensionless scalar field (scalar since charge is an invariant in relativity). There is arbitrariness in the 
definition of e; one can multiply it by a constant and divide all eoi by the same constant without changing anything. 
That is why one must demand that the dynamics of e be invariant under global rescaling of this field (charge-scale 
invariance) . The only possible form of the free action for e is thus 

St = -^lJz-\^(-9) 1/2 d 4 x, (1) 

where g denotes the determinant of the metric g^ u , I is a constant scale of length introduced for dimensional reasons; 
by assumption I cannot be smaller than lp. 

By assumption (1) e must enter into all electromagnetic interaction terms in the matter action via the replacement 
CiA^ i ► eoi eAp, with the usual electromagnetic potential. Gauge invariance of the matter action [assumption (4)] 
will then be preserved only if a gauge transformation means 

eA„ i ^ eA„ + A,„ (2) 

with A any scalar function of spacetime point. In order for the electromagnetic action to be invariant under this 
transformation as well as under e rescaling, it must take the form 

S em = " / F^i-gf/Wx, (3) 

where F^ v = e~ 1 [(eA v ) !l _ l — (eA^)^] obviously stands for the gauge and e-scale invariant electromagnetic field tensor. 
One can consider adding to the integrand of S em a term such as F^F^, where *F flu denotes the dual of F^ v . In 
Maxwellian electrodynamics such an addition is equivalent to a boundary term, and classically irrelevant. Here this 
is not true because of the appearance of the factor e~ 2 in the integrand, but the extra term must nevertheless be 
rejected because it violates time-reversal invariance [assumption (6)]. 



B. Actions 



The appearance of e in the electromagnetic interaction means that the equation for e will involve A^. This was 
found to be a bit inconvenient in Ref. [9]. Thus I shall here adopt the Sandvik-Barrow-Magueijo (SBM) procedure [14] 
of replacing A^ by another 4-potential, a M = eA^. The gauge transformation (2) now turns into 

a M i-» a M + \p, (4) 

so that it is suitable to think of 

ffiiy — Q'l'iix ^^-M (5) 

as the new electromagnetic field tensor; it will turn out to be the physical field tensor (Sec. II D). Likewise I express 
e everywhere in terms of SBM's field ip = lne. 

The total action thus becomes S = S^+S^+S^+Sg, where S m is the matter action (including the electromagnetic 
interaction) 

^(-9) 1/2 d 4 x (6) 

'^U»F v (-9) 1/2 d 4 x (7) 

-g) 1/2 d 4 x (8) 

with 

K = ' „ 8.11 x 10- 26 -^ (9) 
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In the new form of S, tp enters not only in 5*(j) but also appears in the particle masses in S m , even for elementary 
particles. This point is clear from the example of a fermion particle coupled to electromagnetism. If charges were 
truly constant, the process of renormalization would introduce dependence in the fermion (dressed) mass. After 
the transition 

eiAp i ► e 0i eA^ (10) 

the mass becomes a function of e. Swapping for eliminates e from the electromagnetic interaction, but leaves 
e or V dependence in the mass. It is intractable to calculate such dependence by summing the quantum corrections 
to all orders in perturbation theory. However, the framework offers a self-consistent way to compute m(ip) which is 
investigated in Appendix C. 

C. Relation to other theories 

As discussed in Rcf. [9] and already foreseen by Dicke (Ref. [6], Appendix 4), the second form of the framework 
describes constant charges in the presence of varying vacuum permittivity and permeability, e~ 2 ^ and e 2 ^, respectively. 
The framework's action differs from that in Dicke's theory only in that here the electromagnetic Lagrangian is coupled 
to an exponential of the scalar field (as required by charge-scale invariancc of the original form of the framework) 
whereas Dicke made the coupling linear in the field. There is a resemblance also to Jordan's theory of varying natural 
constants [2, 8]; apart from the coefficients of the various terms, the action recapitulated here is the case n — of 
Jordan's. 

This is the place to mention the variable speed-of-light (VSL) theories [20] which have also been touted as variable 
a theories [13, 21]. Dicke's dictum [6] (c.f. Sec. I) that only variation of a dimensionless constant is operationally 
meaningful does not preclude the formulation of a theory which promotes e, or alternatively c, to the status of a 
dynamical field. It only asserts that when e and c appear in the same physical context, the two forms of the theory 
would be experimentally indistinguishable [22, 23]. This point has led to loud controversy with some authors affirming 
that a VSL theory makes different predictions for a variability than a variable e theory [13, 21, 24], while others deny 
it [23]. 

This disagreement is easily defused; c appears in physical actions in at least four contexts: in the electromagnetic 
field-to-matter coupling, cf. Eq. (11), in mass terms, in Lagrangian prefactors, e.g. the he prefacing the Dirac field 
action, and in the Einstein-Hubert action (8). Only for the first of these is c variation fully swappable for e's in the 
context of variable a. In fact, the variable factor in the c featuring in the electromagnetic field-to-matter coupling could 
be absorbed into the electromagnetic potential at the sole cost of introducing a dynamical factor in the electromagnetic 
field action (3), which would then take a form reminiscent of Sf, Eq. (7) of the framework studied here (a variable 
e theory). However, a theory where at least two c's with different roles are promoted to different dynamical fields 
can obviously make predictions different from those of a theory where only e varies. Such, for example, is Magueijo's 
covariant VLS [21] . In the case that all matter is just electromagnetic field, this theory can be understood as one where 
the c of electromagnetism and that from the gravitational action are promoted to different powers of the scalar field. 
This theory which essentially coincides with the general case of Jordan's [2, 8], would predict a different cosmological 
a evolution than does the present framework [13]. But that is because the modified Maxwellian electrodynamics in 
Magueijo's theory is supplemented by a Brans-Dicke style modification of gravity: it is a theory of variable a and 
variable gravitational coupling. 

D. Equations of motion 

Just as in other theories, in the framework the simplest S m is that describing a pointlike classical, possibly charged, 
particle: 

c 7 - 1 (5 3 [x-z(T)]d 4 a; (11) 

where z^(t) = {z°(t), z(t)} is the world line of the particle as function of proper time r, and 7 = dz°/dr is its Lorentz 
factor. The vanishing of SSm/Sx* 1 , when combined with the condition g^dz^dz" = —c 2 dr 2 , gives the equation of 
motion 



mc 



9^- 



az>" az :-\'/' e az^ 



d.r d.r 



(12) 
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or 

Du a dmc 2 e 3 . . 

m-^— = ^(f„ + U a U P 1p t3 ) + —fa(3U P (13) 

where D/dr stands for the covariant derivative along the velocity u a = dz a /dr, itself subject to g^u^u" = — c 2 . It 
is plain from this that f a g is the physical electromagnetic field. Eq. (13) makes it clear that in general, in addition to 
the (suitably modified) Coulomb and Lorentz forces, an anomalous force coming from the tp dependence of mass acts 
on any test charge immersed in a background electromagnetic field. 

The vanishing of 5(S ul + Sif\)/5a^ gives the electromagnetic field equations 

(e-^fn-,u = — f (14) 
c 

f = (15) 

The appearance of e~ 2 ^ in Eq. (14) confirms the interpretation of this factor as permittivity (or reciprocal of the 
permeability) of the vacuum in the present version of the theory. The conservation of the current follows directly 
from Eq. (14); the conserved charge is the truly constant e . It is in this sense that electric charge is still conserved 
in this "variable charge framework" [9]. 

Finally the vanishing of S(Sm + S^) + S(f))/6ip gives the equation for ip, 

W = 4^^^i^l - £e-* W - (16) 
dip IV ~9 2 

while that of S(S m + + S(f) + S g )/Sg^ v gives the gravitational field equations. 

To fully specify the equation of motion of charged particles, Eq. (13), one must specify m(ip). As mentioned 
in Sec. I, it has been customary to assume that the spacetime dependence of mass reflects the electromagnetic 
contribution to it, here proportional to e 2 or e 2 ^ . Thus for ordinary matter (for which the Coulomb energy far 
surpasses magnetic energy), it has been customarily assumed that a fixed fraction ( of mc 2 is Coulombic, so that 
dm/ dip = 2(mc 2 [9, 10, 11, 12, 13, 14, 15, 16, 17, 18]. Because C is expected to vary from object to object, violation 
of the WEP would seem to be inevitable. However — and this is one of the main points of the present paper - 
the mentioned prescription for m(ip) is the wrong one for the framework defined in Sec. II. The point is that the 
field equations (14)-(16) are nonlinear, e.g. the permittivity e -2 ^ is determined by the electromagnetic field strength 
through Eq. (16). It is well known that the equations of motion of the sources of nonlinear field equations, e.g. 
Einstein's equations, cannot be freely prescribed. Rather, their nature is specified by the field equations themselves. 
Here, too, one must let the field equations specify the nature of m(ip) which defines the explicit form of the source's 
equation of motion, Eq. (13). 

III. ISOLATED POINT ELECTRIC CHARGE 
A. The charge's fields 

An important step in clearing up the status of the WEP in the framework is the study of the motion of a pointlikc 
charge in a specified external electric field. The theory's nonlinearity does not permit one to ignore the charge's own 
fields in setting up its equation of motion; the task would be easier in Maxwellian electrodynamics where fields can 
be superposed. For now I neglect the curvature of spacetime; as shown in Sec. HID, this is entirely justified for 
elementary particles and small collections of them. I also assume that initially the charge is held in place by some 
unspecified force, and then released, so that its fields are static to start with. 

Setting 7 = 1 and x(t) = in Eqs. (14)-(16) one can look for a static solution with vanishing magnetic field 
{f ij = 0). With the notation c^E = {/ 01 , f 2 ,.f 03 } one gets 

V- (e^E) = 47re <5 3 (x) (17) 
,rdmc 2 „ q , . 1 



V 2 ip = 4ttk 2 



6* ^) + I^ E 1 (18) 



By Eq. (5) and the obvious condition that a M may be taken time independent, E must be a gradient: E = VT. 
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The generic solution of Eq. (17) is 

e -2V V T = -V$ + b; V-b = 0, (19) 

where 

$ = e /r (20) 

is the usual Coulomb potential of the charge eo in spherical polar coordinates {r, ip}. A particular solution is 
obtained when b equals some constant vector <E. Since ^ should asymptote to a constant at infinity, (£ is obviously 
the applied external electric field up to a positive proportionality constant. The curl of Eq. (19) gives x VT = 
which shows that ip and T are functions of each other. By Eq. (19) both can be taken as functions of only the potential 
V = $ - <£ • x. Thus 

E = -e^'VV (21) 

with tp = ^>(V). So far it is clear that Eq. (21) is the unique solution when b = <£. But as argued in Sec. IV A, the 
uniqueness survives when b is only known to asymptote to (E. 
In view of Eq. (21), 

V 2 ^ = ^'V 2 y + V"e~ 4,/J E 2 (22) 

(' = d/dV, etc). Since V 2 V = — 47reo(5 3 (x), the first (second) term on the r.h.s. of Eq. (22) must match the first 
(second) term on the r.h.s. of Eq. (18). The second condition, namely, 

V>" = n 2 e 2 ^ (23) 

can be integrated after multiplication by ip' to give 

V/ 2 = K 2 (e 2 ^+vu) (24) 

with vo a dimcnsionless constant of integration. Appendix A obtains all solutions of this equation. 
Identifying the pointlike source terms of Eqs. (18) and (22) gives 

K 2 (dmc 2 /dip) = ~e tp' at x = 0. (25) 

This condition makes the two source terms equivalent even when the 5 function is somewhat smeared (I shall show 
below that in this theory there is a lower bound to the radius of any charge, so this smearing must occur). As evident 
from the appearance in the r.h.s. of a term linear in eo, the dependence arising from this equation is not automatically 
of the form m = m + const, x e 2 e 2 ^ commonly adopted (Coulomb model). 



B. Physical choice of vo and WEP 

To find out about w I now work out the total force of electric origin acting on the charge eo. I make no attempt 
to separate out the self field. Immediately after the charge is released Eqs. (13) and (21) give 

mdv/dt = -(dmc 2 /dip)Vip + e 'E 

= -[(dmc 2 /diP)i>' + e e 2 ^]VV. (26) 

It is to be stressed that the force on the r.h.s. here comprises both the (modified) Coulomb force and the anomalous 
force. Substituting (dmc 2 /dip) from Eq. (25) and ip' 2 from Eq. (24) gives, after a cancellation, that 

mdv/dt = zue W (27) 

The choice w = in the solution corresponding to a charge eo is physically untenable. It would mean that in 
the presence of an arbitrary external field <£, eo experiences no force whatsoever. This is contrary to all experience, 
and the theory developed here is supposed to describe the real world. It must be, then, that w ^ 0. For vo ^ 
the force on the r.h.s. includes the usual self-force term proportional to V$ which is always dropped in Maxwellian 
electrodynamics; I do so too. Thus 



mdv/dt = —zue (£ 



(28) 
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It would be hasty to reject outright all the cases with zu 7^ — 1 on the grounds that Eq. (28) would then have the 
wrong form; after all <£ is the external electric field only up to a positive proportionality constant. However, with 
zu > the net force acting on the particle would be opposite the accepted one. One cannot remove this problem by 
assuming that the "passive" charge here which senses the applied field is opposite in sign to the "active" charge eo 
which is the source of the electric field in Eq. (17). For if this were true here, it would be true in other cases too. 
If (£ is then interpreted as representing the field of a distant charge (up to a positive multiplier), it is immediately 
apparent that charges of one sign would attract. I thus conclude that physically zu must be negative. 

In Appendix A all solutions of Eq. (24) are found; that for negative zu is (x is an integration constant): 

e^ = M 1/2 scc(\zu\ 1/2 kV + X ) (29) 

This solution is unique (in the physical sense) by the one-to-one correspondence between solutions of Eq. (24) for one 
charge and the solutions of the corresponding equation for many charges (see Sec. IV A). The latter are certified as a 
complete set of solutions by the duality argument to be set forth in Sec. VIA, and the particular solution with zu < 
corresponds to Eq. (29). 

Let us interpret e p = Im^^eo as the physical charge and <£. p = \zu\ x / 2 <£ as the physical external field. This brings 
Eq. (28) to precisely the everyday form of the Newtonian equation of motion in an external field; zu disappears from 
the equation of motion. 

Eq. (29) now takes the form 

= \vj\ 1/2 sec(KV p + x)- (30) 

where V p = $ p — <£ p • x, and <£> p is built just as $ in Eq. (20) but from the physical charge. It would seem that zu still 
appears after the reinterpretation. However, this is just an illusion. Consider the energy density of the electric field, 
£ = (87r) _1 e _2 ^E 2 according to Eq. (B2). This can be recast in terms of physical quantities in such a way that zu 
does not appear explicitly: 

£ = (Sir)" 1 sec 2 ( K V p + X ) ■ (« P - V$ p ) 2 (31) 

The same can be said about the purely Coulomb force on the test charge, $ = e E, which takes the form 

$ = e p sec 2 (nV p + x)-(<£ P - V$ p ). (32) 

Again zu has disappeared. All this means that setting zu = — 1 (identifying eo with the physical charge) is a matter 
of convention, and does not entail any assumptions about the physics. I thus take zu — — 1 and drop the subscript p 
henceforth. 

The equation of motion (28) for the particle in an external field (with zu = — 1) already shows that the WEP is 
not violated — at least not in an obvious way — in the theory in question. The putative anomalous force from the ip 
dependence of the mass is compensated for by the modification of the Coulomb force effected by the ip ne ld [the e 2 ^ 
factor in Eq. (26)] in such a way that, at least for a quasistationary charged particle, the force acting on the charge 
is independent of the fraction of its rest mass which is of Coulomb provenance. One finds no indication here for a 
violation of the WEP. 

A brief remark about neutral particles is in order here. In that case Eqs. (17)-(18) reduce to a Poisson equation for 
tp with a pointlikc source. The formal solution is ip — C/r with C a constant proportional to dm/ dip of the source. 
Dvali and Zaldarriaga [12, 25] suggest that C is generically nonzero with the consequence that long-range anomalous 
forces exist even between neutral particles. For the neutron they infer this from the a dependence of the neutron 
mass coming from virtual photon exchange between its constituent quarks [12] as calculated perturbatively to low 
order. As discussed in Sec. IV A, the framework suggests rather that C — 0, at least when some charged particles are 
also present. 

C. Minimal size of an isolated charge 

Henceforth I set <£ = 0. With zu = -1, Eq. (30) gives 

= sec(K$ + x)- (33) 

Apart from the inclusion of the "phase" x> this is the original solution presented in Ref. [9], where it was pretty much 
gotten by guessing. Redefining \ is equivalent to shifting the zero of $; this is certainly without physical significance 
here as in Maxwellian electrodynamics. As already clear from Eq. (20), I here adhere to the convention that <i> 
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vanishes at infinity. The value of \ is thus fixed by the asymptotic value of , which coincides with the instantaneous 
cosmological value of in the appropriate model of the universe. 

According to Sec. II, it is permissible to multiply all eoi by a common positive constant while simultaneously 
dividing e by it. I exploit this freedom to set the cosmological value of e to unity at the present epoch. After this is 
done one can define -0 as the logarithm of e and pass to the second form of the theory. From the solution (33) at the 
mentioned epoch one then finds that \ = (this was also the choice of Ref. [9]). Although x evolves cosmologically, 
it is evidently possible to set x = at any one cosmological epoch. Thus results we shall obtain below that depend 
on having \ = are valid at any single epoch. It is only when interests centers on comparing physics at two separate 
epochs that one can no longer do away with \. 

By Eq. (33) can diverge and then turn negative when |$| > 7r 3 / 2 (flc) 1 / 2 ^ - 1 . Because the permittivity of the 
vacuum cannot be negative, this must mean that |$| can never reach such values, i.e., that the particle with charge 
e is spread over a sufficiently large radius R to prevent this. By Eq. (20) this condition is 

R>TT- 3 / 2 (eo 2 /hc) 1 ^(l/£p)e P . (34) 

Because asymptotically — > 1 at the present epoch, one can take (eo 2 /he) 1 / 2 = (137) -1 / 2 <~ 0.1 for e the 
elementary charge. And by the framework's assumption (7), I > £p; in fact I shall show in Sec. VII that if the alleged 
cosmological a variability is to find explanation in this framework, I must be an order of magnitude above £p (see 
also Refs. [9, 13, 14, 15]). Thus the lower bound on the radius of any charge is at least a Planck length. Composite 
particles, e.g. the proton, easily satisfy Eq. (34). For leptons and quarks which are regarded as pointlike, quantum 
gravitational effects must intervene at radii of a few Planck lengths and modify the above classical considerations. But 
it is noteworthy that our formal lower bound on R is not at variance with the widespread belief that no elementary 
particle can be smaller than £p, the scale at which spacetime can no longer be regarded as a continuum. 



D. Why neglect spacetime curvature ? 

The neglect of spacetime curvature in all the preceding calculations may be justified when the particle in question 
is either elementary or made up of a not excessive number of elementary particles. Appendix B shows that corrections 
to the usual exterior Reissner-Nordstrom metric (of general relativity) belonging to a particle with charge eo become 
important in our theory only when |$| is no longer small compared to ^^(ftc) 1 / 2 / -1 . By Eq. (34) |$| gets that big 
only as r approaches R. Thus, if all one wants is to investigate the source's exterior, one can, with good accuracy, 
employ the Reissner-Nordstrom metric, 

ds 2 = -e- x dt 2 + e x dr 2 + r 2 (d9 2 + sin 2 6dy 2 ) (35) 
e- x =l-2Gmc 2 r- 1 +Ge a 2 c- 4 r- 2 , (36) 

where m is the source's mass. 

This metric begins to departs seriously (5% level) from Minkowski's at a radius r <~ max(5G 1 / 2 eoc -2 , 40Gtoc~ 2 ). 
Now for the known charged elementary particles or small agglomerates of them, m <C G _1 / 2 e (in fact G _1 / 2 eo/m is 
of order 10 25 for the electron and 10 20 for a nucleus or atom). Thus significant departures from flatness are only found 
at r < 5G 1//2 e c~ 2 = 5(eo 2 /hc) 1 / 2 £p, i.e., at Planck scale where the whole classical description is already irrelevant. 
And as pointed out in Sec. IIIC, pointlike particles in this framework cannot be smaller than this. Thus, description 
of the exterior of a pointlike elementary object (or a small collection of such) can well afford to ignore spacetime 
curvature. 



IV. MULTIPLE ELECTRIC CHARGES 



A. The solution 



The results above may be generalized to a collection of N charges eoi initially clamped at positions Zj, i = 1, 2, • • • , N. 
Again with neglect of gravity Eqs. (14-(16) reduce to 



VV = 47TK 2 



V- (e" 2 *E) =47r^e 04( 5 3 (x-z l 

i 

E^ 3 (— ) + ^E 2 



(37) 
(38) 
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As in Sec. Ill, I take E = VT. 



I define anew 



*(*) = £ 



X - Zj| 



(39) 



this is the standard Coulomb potential due to all the charges. Recalling that 



V 2 $ = -47r^e M (5 3 (x-z 4 ) 



(40) 



it is now easy to check that 



E = -e 2,/, V$ 



(41) 



together with any of the choices (x is a constant) 



( ±(k$ + x)~ 1 ; ^ = 

e^ = < ± ^/n7CSch(^/n7K$ + x); w > 
I y/ \zu\ sec(y/ \zu\ k& + x); tn < 



(42) 



constitute solutions of Eqs. (37)- (38) provided the analog of Eq. (25) is satisfied for every one of the charges, to wit 



I shall exploit this last equation in Appendix C to elucidate the spatial variation of the m, . 

It turns out that Eqs. (41) and (42) comprise all the static solutions for multiple charges. Although it is hard to 
prove this directly, results in Sec. VIA enable this be established immediately by duality arguments. The one-to-one 
correspondence between the branches of ip in Eq. (42) and those for the one-particle solution (A5) then establishes 
that in Sec. IIIB, Eq. (21) is indeed the most general form for E possible even if b is not assumed constant except 
asymptotically; generalizations of Eq. (21) would necessarily affect tp thorough Eq. (18). 

Since the physical multiparticle solution should include the physical single particle solution as a special case, I 
put w = — 1 here as in Sec. IIIB. I also set x = on the basis of the argument of Sec. IIIC. That this choice of 
parameters is consistent is made clear by the following argument. Suppose $j>i, the part of <E> coming from particles 
i = 2, 3, • • • , N, is approximated by the first two terms of a Taylor series about zi. Denoting «$i>i(zi) by xi an d 
V<f>i>i at zi by — <£, the last of Eqs. (42) gives in the vicinity of charge i = 1 



This is of the same form as Eq. (30) for the single charge solution. Although the phase x of the whole charge complex 
has been set to zero, each individual charge has an associated phase, Xii induced by its neighbors. The result (44) 
makes it clear again that the one-charge solution Eq. (30) is physically unique, as mentioned earlier. 
With the above choices the full multiparticle solution is 



This result serves to clear up the question (Sec. IIIB) about the value of dm/ dip for a neutral particle. Because 
neutrinos and neutral mesons are scarce in laboratory matter, the particle of most interest is the neutron. It is 
composite, but I avoid any discussion of its structure and extension, and treat it as pointlikc object as it would indeed 
appear to low energy probes. In this spirit one can include a neutron as particle n in the collection just discussed by 
formally taking eo n — » 0. It is easy to verify that this limit is a solution of the equations. It follows from Eqs. (45)-(46) 
that tp is regular at the neutral's position, and from Eq. (43) that dm n /dtp — 0. More generally, for a neutral pointlike 
particle, m is tp (and spacetime) independent. It remains a task for the future to reconcile this conclusion with the 
dependence of the neutron mass on a according to perturbative calculations within QCD [12]. 



k 2 (dm-iC 2 1 dtp) = —eoiip' at x = z,. 



(43) 




(44) 



e7 = sec(«<I>) 

E = -e 2 ^V<5> = - K - 1 Vtan(/«$) 



(45) 
(46) 



B. The momentum equation and WEP 



In our context the WEP would require that a composite possibly charged body moves in uniform external electric 
and gravitational fields with an acceleration which depends only on its total mass and charge, but not on its detailed 
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inner structure. To put the framework to the test in this respect, I imagine all N charges in "the world" to be lumped 
into two clusters. One, a spherical massive one, I put at the origin and regard it as a single particle of mass mi and 
charge e i- The other, denoted by C, is made up of N — 1 charges eoi with masses moi; i = 2, 3, . . . N. By assuming 
the distance between the clusters is large compared to both their extensions, it is possible to think of the charge at 
the origin as pointlikc, and C as immersed in the uniform fields of the former. 

I take the static solution for the fields, Eqs. (45)-(46), as part of the initial conditions for the envisaged dynamical 
situation. I restrict discussion to the initial moment t — when all the charges are still at rest. This restriction is 
necessary here because the fully dynamical solution is as yet unknown. However, it should be clear that any violations 
of the WEP would be expected to show up already at the initial moment because they involve the acceleration. The 
demonstration in Sec. IV C that no such violations occur, at least within the approximation to be described presently, 
strongly suggests that the WEP holds to great accuracy in the dynamical situation as well. 

As a first step I calculate the rate of change of total momentum in the masses constituting C. According to Eq. (12) 
in the nonrclativistic approximation and with neglect of gravity 

J N N Q 2 

5>v 4 = £(- W + e „,E) . (47) 

»=2 i=2 r 

Substitution from Eqs. (43) and (45)-(46) and use of the identity sec 2 x — tan 2 x = 1 transforms this into 

N N 

(II 



d TO * V * = -J2 e o* v *lx= Zl (48) 



Now near charge i, V$ is dominated by the self-field (x — z»)/|x — Zi\ 3 which points radially out from it. As in 
Maxwellian electrodynamics, here one may think of the force from this part as averaging out in the limit. Left over 
from the force are the terms (j ^ i) 

y y e oi eoM - z 3 ) = y e^z, - Zl ) 

with the second form following cancellation of the j ^ 1 terms in pairs. The assumed smallness of C as compared 
with the distances |z, — zjj justifies replacement of every Zj on the r.h.s. by the cluster's center-of-mass position Z 
(an approximation no different from the one customarily made in Maxwellian systems): 

^ m ' V ' = Q |Z- Zl| 3 ( 50 ) 

with Q = X)il=2 e 0i- Thus the rate of change of momentum of the particles in the cluster is controlled by the 
formal Coulomb field of charge i = 1, approximated as uniform at C, and by C's total charge. Just as in Maxwellian 
electrodynamics, here the cluster's internal structure does not affect the rate of change of its total particle momentum. 

Although this finding is consistent with the claim that the WEP is satisfied here, it is no complete proof: the 
rrii vary, making the relation between rate of change of the total particle momentum and the acceleration of the 
center-of-mass' less clear than usual. 

C. Restricted microscopic proof of WEP 

To clarify the above point I look at the time component of energy-momentum conservation for the whole system (C 
and charge i = 1). If T^ v includes the tensors for the particles and the electromagnetic and scalar fields, T^ ; „ = 
can be derived from the gravitational field equations as usual. Ignoring the gravitational field gives 

dT oa /dt + dT 0l /dx* = 0. (51) 

Multiplying the equation by x- 7 , integrating over a large volume V containing C but excluding charge 1, and integrating 
by parts gives 

^ / T°Vd 3 x = / T°3d 3 x - I xiT 0l d 2 S t . (52) 
dt J v J v Jq V 
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The integral on the l.h.s. is just the product of total mass M = J v T 00 d 3 x of the cluster (neglecting contributions 
from the fields beyond the reaches of C) and the component of Z. One further time derivative gives 

^ / T 00 xid 3 x=4- I T°^d 3 x-^I x^T 0i d 2 Si. (53) 
dt 2 J v dt J v dt J av 

The spatio-temporal component T° l figuring in Eq. (53) receives three contributions. The first, T)/) 01 from the 
field (see Appendix B), comprises products of f° l with pK The evidently vanish at t = (no motion, no magnetic 
field). Further, because I assume the static solution (45)-(46) for ip and E holds at t = 0, foi = —cE l is a gradient. 
Then the identity + foij + fjo,i = (f^ v derives from a 4-potential) shows that initially = 0. Therefore, time 
differentiation anihilates both integrals over Tm ' in Eq. (53) at t = (that ip = at t = will be shown presently). 

The second contribution, T^ 01 (again sec Appendix B), comprises products of V' 4 with ip. Now because the 
scalar equation (16) is of second order in time, one may require ip = at t = 0; this is consistent with this instant 
representing the end of a purely static situation. Further, comparison of Eqs. (38) and (16), the latter satisfied 
identically by solution (45)-(46), shows that also ip = at t = for the assumed initial conditions. Thus time 
differentiation of both integrals over T^° l in Eq. (53) fails to produce nonvanishing contributions at t = 0. 

Accordingly, the time derivatives at t = in Eq. (53) of both integrals over T°* come solely from the particles in C. 
Since the surface dV lies beyond them, the r.h.s. of Eq. (53) comprises solely J T m 0j 'd 3 x. The tensor in question is 
derived from action (11) by variation of g^ v and use of g^u^u" = —c 2 : 

Neglecting the difference between r and t, replacing y/—g by its Minkowski value c, and substituting in Eq. (53) gives 

dt 2 - dt 2^ m ^- ^ 

I complete the derivation by showing that at t — M can be taken out from under the derivatives. Integrating 
Eq. (51) over V and using Gauss' theorem gives 



dt J v 



T d x = -<p T^d z Si. (56) 

«9V 



Differentiating this result by t gives 



dt 2 



[ T m d 3 x = ~<f T 0l d 2 S t . (57) 

JV dt Jg v 



By precisely the same arguments given above, the r.h.s. of these two equations vanish at t = 0. Since the integral in 
both l.h.s. is Cs total mass M, this demonstrates that M = M = 0ati = 0. Of course this does not mean that M 
is a conserved quantity, only that it behaves as M(t = 0) +0(t 3 ) for short times. Substituting this result in Eq. (55) 
and taking Eq. (50) into account gives at t = 

~ _ Q egi(Z-£i) 

Z -M |Z- Z1 |3 ' (58) 

Thus in harmony with the WEP, the acceleration of Cs center-of-mass in the field of the distant charge eoi is fully 
determined by its mass M and total charge Q, and is insensitive to its structure (disposition of the member charges, 
their charge to mass ratios, etc.) To reach this result the mass M of C had to be identified with the integral of T 00 
taken over a finite, albeit large, region. This is a necessity in any situation when the system of interest is not the 
only one in the universe; the same procedure would be required in any other field theory. Another limitation of the 
approach is that result (58) is rigorously valid only when all particles in C are assumed to be at rest (thus the approach 
neglects purely magnetic effects). But intuitively the acceleration's universality property should remain valid if all 
velocities are small and spin magnetism is weak, as I indeed show in Sees. VB and VI B. 

The results here concur with those reached by Landau, Sisterna and Vucetich [19] on the basis of the THe/z formalism, 
a nonrelativistic generic parametrized microscopic description of the gravitational, particle and electromagnetic sectors 
of field theories. The authors find a connection between the lack of overall charge conservation and violations of the 
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WEP. Looking at the present framework in the light of experimental constraints on charge nonconservation, they are 
able to certify that WEP is respected to a fractional accuracy (in the acceleration) <~ 1CP 18 , beyond the projected 
sensitivity of WEP tests in the foreseeable future. Since charge conservation is actually exact in the framework 
[see the comments after Eq. (15)] the implication is actually stronger: no WEP violations are expected, at least 
nonrelativistically, even when slow motion of charges is allowed. 

There exist other possibilities for variability of a by way of a scalar field which do not run counter to the EDB 
experiments. Such is the supersymmetric grand unified theory of Chacko, Grojean and Perelstein [26], according to 
which a late epoch cosmological phase transition causes a jump in a while generating a vacuum expectation value 
of the scalar field which makes it short ranged. Consequently, although a charged particle's mass can be scalar field 
dependent, this does not lead to long range anomalous forces which would contradict the said experiments. 

V. MACROSCOPIC DESCRIPTION OF ELECTRIC STRUCTURE AND WEP 

A virtue of the "proof" of the WEP in Sec IV C is that it does not require knowledge of the equations of motion 
for macroscopic matter [as opposed to the microscopic level Eq. (13)]. Its shortcomings are that it establishes the 
validity of the WEP only for brief time intervals after a putative quiescent situation, and only when the magnetic 
structure of the systems involved is ignored. Both shortcomings will now be remedied. By averaging over microscopic 
quantities to produce a macroscopic (or continuum) version of the theory, I show in Sec. VB that the WEP-breaking 
effects pointed out in Refs. [6, 7, 9, 10, 11, 12, 13] are ruled out if one is willing to assume that a macroscopic chunk 
of matter moves according to Eq. (13) with a natural definition of m(ip) obtained by macroscopic averaging. No 
short-times assumption is then necessary. And in Sec. VI B magnetic dipolc structure of matter is incorporated into 
the arguments. 

Sees. Ill andIV describe matter microscopically, that is as a collection of pointlike particles i with definite charges 
eoi (some of which may be zero) and masses nii subject to the relation (43). One key assumption I make on the way to 
the macroscopic description is that the many charges solution (45)-(46) remains valid, at least approximately, when 
the particles move slowly. I define the spatial average Q of a quantity Q as the integral of Q over some macroscopic 
region of volume V, all divided by V. In flat spacetime or in local Lorentz frames (assumed to be large enough to 
encompass the said macroscopic region), no ambiguity in this definition arises from issues of parallel transport of 
vectors and the like. 

Another assumption I make is that for a macroscopic test body T of mass p and charge q moving on background 
scalar, Newtonian and electric fields, jp, V</>jv and <£ (all regarded as approximately uniform), moves according to 
Eq. (13) [Eq. (26) in the nonrelativistic case]: 

pdw/dt = -pV<f) N - {dpc 2 /d^)V^p + q<£ (59) 

The assumption would be a triviality but for the stipulation that p(-tp) is to be identified with the macroscopic average 
of the total energy density in T multiplied by its volume. Evidently q is to be interpreted as the sum of charges in 
T, while (£ is the electric field determined from its sources by Eq. (37) and averaged over the volume of T. The only 
quantity which requires special discussion is V-0. 

A. The Coulomb model for «/> 

In the language of the present framework, Dickc's argument [6] amounts to assuming that the macroscopic field tp 
is related to its source S via a macroscopic version of Eq. (38) with the sum of (drriiC 2 / dip)5 3 (x — z^) replaced by 
2(pc 2 , where £ is the typical fraction of the source's mass density p which is of Coulomb provenance. This comes from 
assuming that pa a a e 2 ^. Now the second term in the square brackets in Eq. (38) is, according to Appendix B, 
exactly twice the density of electric energy, just as the former contribution. Hence in the Coulomb model for the 
source, Eq. (38) takes the form 

V 2 ^ = 16ttk 2 Cpc 2 = (l/£ P ) 2 ((/nc 2 )4TrGp, (60) 

A likeness of this equation and/or its cosmological version reappears in subsequent treatments [9, 10, 11, 13, 14, 15, 
16, 17, 18]. 

The similarity of Eq. (60) with Poisson's equation for the Newtonian potential of <S, and the similarity between the 
asymptotic boundary conditions on 4*n and ip, permits the identification 



V0 = (l/lp) 2 ((/nc 2 )V<f> N 



(61) 
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Actually, in a nonspherical configuration the curl of some vector may be added to the r.h.s.; however, this "correction" 
must decay asymptotically as 1/r 3 (cf. Ref. 27, Appendix A) and so should become irrelevant some distance outside 
S. Eq. (61) shows that Va/« = (1/£p) 2 (2(/ttc 2 ) V</>at. With I ~ ip this result is very much like Dicke's conjectured 
one [6] mentioned in Sec. I, except for the replacement of a by 2£, which for ordinary matter is of the same order as 
a. 

How does this all bear on the WEP ? Suppose T moves in the vicinity of S. In the Coulomb model one replaces 
(dfic 2 /dip) in the equation of motion for T, Eq. (59), by 2('/j,c 2 , where (' is now the Coulomb energy fraction of jic 2 . 
With the replacement (61) this gives (in the absence of other sources) 

dv/dt = -[1 + (l/£ P ) 2 (2(C /ti)]V4>n + q£, (62) 

a version of which first appeared in Ref. [9]. Now, as mentioned in Sec. I, for ordinary matter £ is of order a few 
times 10~ 3 and varies by about 10~ 3 from material to material. The latest EDB tests of WEP find that dv/dt is £' 
independent to fractional accuracy 10 -13 [8]. This is consistent with Eq. (62) only if I is 10 _3 ^p or smaller. Because 
of this and assumption (7) of the framework, I inclined in Ref. [9] to the opinion that there is no a variability in 
nature (as would be the case if I = 0, e.g. electrodynamics exactly Maxwellian) . 

This was also the conclusion of Livio and Stiavelli [10] who noted the difficulty in explaining the alleged cosmological 
a variability with I as small as 10~ 3 £p and the accepted matter content of the universe. Olive and Pospelov [11] also 
took cognizance of this problem; to solve it they proposed that the cosmological dark matter is much more strongly 
coupled to the ip field than is ordinary matter. This would have the effect of counteracting the smallness of l/ip 
inferred from tests of the WEP. By contrast, Magueijo, Barrow and Sandvik [13] see no immediate strong contradiction 
between the EDB experiments tests; they infer {l/ip) 2 w 10~ 4 from the claimed cosmological a variability by assuming 
that cosmological matter has ( w 1. They further adopt, for ordinary matter, the very low values ( <~ (' ~ 10~ 4 . 
However, these run on the face of simple estimates from nuclear Coulomb energy (see Sec. VB2). There is thus a 
tension between the claimed cosmological a variability and the tests of the WEP. In the next subsections I show that 
the blame for the above impasse lies squarely with the misuse of the Coulomb model, there being no need to take 
I < ip. When needed in what follows I assume I is above ip and within an order of magnitude of it. 

B. Coulomb energy and WEP violation 

For a quick orientation one can describe the detailed electric structure of the source S of ip with the static electric 
solution, Eqs. (39), (45-(46). One is thus neglecting effects of internal motions and magnetic structure, both usually 
minor complications for nonrelativistic sources of interest. 

I assume that <C 1 everywhere and check this on a case-by-case basis. Then Taylor expanding both sides of 
Eq. (45) and averaging as stipulated in the preamble of this section gives 

$ „ I K 2^2 = I K 2^2 + l K 2/ $ _ ^2 (gg) 
2 2 2 

Note that matter with electric structure causes the physical a, which is proportional to e 2 ^ « 1 + 2tp, to be slightly 
larger nearby than asymptotically (using the Coulomb model Refs. 13, 14 predict an effect with the opposite sign and 
much larger magnitude). Now outside a macroscopic source S the fluctuation term here should be relatively small 
compared to $ 2 , except in rather artificial situations where $ very nearly vanishes (the net charge and some higher 
multipoles are exactly zero). The latter are not important in our context, so henceforth I shall drop the fluctuation 
and write just $ for $. 

1. Natural sources of gravity and ip 

First I look at natural sources of e.g. the Sun and Earth in the EDB experiments. Unless S is (almost) exactly 
neutral, $ is dominated by its monopole part. For example, Earth is known to bear a net charge at any time. Thus 
for a quasispherical source $ w Qr^ 1 and |V^| w n 2 Q 2 r~ 3 at distance r — |r| from the source's center. 

There is a natural bound on Q if the source, a natural large object like Sun or Earth, is to be quiescent: \Q\ < 
GMrrip/eop, where M is 5's mass, and m p and eo P the proton's mass and charge, respectively. For if Q > GMm p /eo P , 
the source's electric field can drive away any free protons formed nearby by, say, cosmic ray ionization of hydrogen 
even against the pull of gravity, while S captures the electron and so decreases Q. And if Q < —GMm p /eo Pl S 
can certainly drive away the free electrons, capture the protons and so decrease (In the above argument it is 
important that e 2 ^ w 1 by assumption, soEk — V<&.) 
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In view of the restriction on Q, 



l^l / 1 M/2/ 1 slpGM 



°Q P /hc w 0.0073 and A p = h/m p c w2x 10~ 14 cm is the proton's Compton length. Obviously GM/c 2 r < 1 
because the source is not a black hole. It follows that < 3 x 10 (l/£p) <C 1, thus verifying the initial assumption. 
Further, GM/r 2 w |V^>jv|- Thus |V^| is a fraction 2.2 x 10~ 37 C _1 or smaller of the Coulomb model prediction (61). 
In view of Eq. (62), violations of the WEP reflecting S's electric structure show up only at the acceleration fractional 
level 10~ 40 (l/£p) 2 or smaller. (In light of the conclusions of Sees. IIIB and IV C, the tiny residual WEP violation 
here, and similar ones below, may well be artifacts of the macroscopic averaging procedure.) As mentioned already 
the extant laboratory EDB experiments are sensitive only to acceleration fractional differences at the level of 10~ 13 . 
The planned satellite borne STEP experiments are projected to be sensitive only at the level 10~ 18 . 

What if S is nearly neutral, so that $ is dominated by its dipole moment p: $ w p • r r~ 3 . In that case |p| should 
be on the order of the radius R of S times a typical separated charge \Q S \- The magnitudes of the charges ±Q S , which 
pictorially reside on opposite polar caps of S, are restricted by the same inequality as the monopole charge above. 
Otherwise free electrons near the negative polar cap would be driven away against gravity and conveyed by the dipole 
field to the positive cap, thus helping to diminish |p|. Repeating the above argument one finds the expressions for 
k|$| and |V^| to be similar to those in Eqs. (64) but each with an extra factor (R/r) 2 < 1. Thus k|$| <C 1 is still 
satisfied, and the failure of the Coulomb model is accentuated. Evidently higher multipoles do not offer a way out of 
the conclusion that the Coulomb model is very far off the mark. Evidently laboratory and space tests are far from 
sensitive to violations of the WEP coming from the electric structure of the sources of gravity. 



2. Laboratory sources of gravity and ip 

But what if the source S is not a gravitating body in the dirty interplanetary environment, but rather a mundane 
body in a clean laboratory where its charge docs not immediately get neutralized ? Would the Coulomb model apply 
then ? An example might be furnished by a lead sphere of radius R. Again Eq. (63) tells us that if the sphere 
holds charge Q, then outside it |Vi/>| w n 2 Q 2 r~ 3 < k 2 (Q 2 / R)r~ 2 . The factor Q 2 /R is twice the formal Coulomb 
energy associated with the net charge Q. How big can Q be here ? Unless the sphere is in an evacuated cell (which 
procedure can buy us a few orders of magnitude in the discouraging results below) , it can only be charged until its 
surface electric field reaches the air breakdown level, ~3x 10 4 V/cm = 1 x 10 2 esu. At that point the sphere is at 
$~lx 10 2 csu R. According to Eq. (9), k w 8.11 x 10~ 26 esu. Thus indeed k|$| <C 1, as assumed. In said state the 
sphere holds ~ lx 10 2 esu R 2 of free charge, and thus twice its macroscopic Coulomb energy is <~ 1 x 10 4 erg/cm 3 R 3 . 
Hence 

\Vtp\ < 1 x 10 4 erg/cm 3 K 2 i? 3 r~ 2 (65) 

By contrast, in the Coulomb model Eq. (60) would predict 

|V$ = 4 K 2 [C(47ri? 3 /3)pc 2 ]r- 2 , (66) 

with the square brackets recognizable as the microscopic level Coulomb energy. Being lead (p = 1.13 x 10 1 g/cm 3 and 
mass number 207.2) the sphere contains ~ 1.38 x 10 23 /cm 3 i? 3 Pb nuclei; each contributes ~ 795 MeV of Coulomb 
energy [29] for a total of <~ 1.75 x 10 20 erg/cm 3 R 3 . Thus in the present example the true |V-0|, Eq. (65), amounts 
to a fraction 1.4 x 10~ 17 of the Coulomb model prediction, independent of R. In this case also the Coulomb model 
yields a resounding overestimate. In view of Eq. (62), violations of the WEP for test body motion subject only to the 
field of a macroscopically charged laboratory sized object are at a fractional level 10~ 23 (l/£p) 2 . But this fraction is 
actually suppressed to 10 _31 (//£p) 2 because S's gravitational field is diluted by a factor 10 8 by Earth's gravitational 
field (for this calculation I take <S's mass to be 10 6 g). Thus the WEP violation considered here is well below the 
sensitivity of the laboratory EDB experiments. The planned STEP experiments are irrelevant in this context because 
they do not include their own source of ip. 

The prospects for WEP violation improve substantially when macroscopically charged sources in the laboratory are 
replaced by macroscopically polarized ones in Earth orbit (an extension of the STEP experiment which is planned to 
carry only test objects). A good example is a sphere of radius R made of ferroelectric material. In a ferroelectic there 
is spontaneous alignment of the molecular electric dipoles, so that a macroscopic polarization vector P appears. Thus 
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in Eq. (63) $ s=s (47ri? 3 /3)P • r r~ 3 . A good estimate of the maximum |P| is one elementary dipole eo p d per molecular 
volume d 3 . With d <~ 10~ 7 cm this works out to |P| < 4.8 x 10 4 csu. Indeed PbTi03, a ferroelectric with one of the 
largest measured polarizations, shows |P| ~ 1.5 x 10 5 esu [30]. Given that r > R, one finds < {Att)^ 1 / 2 a 1 / 2 lRd~' 1 . 
With R <~ 10 2 cm this bound is below 10 -18 (l/£p), small as originally assumed. 
Using Eq. (63) one estimates 

|V0| ~ 4 K 2 (i?/r) 3 [(47ri? 3 /3)P 2 ]r- 2 (67) 

Eq. (67) exhibits an energy P 2 rf 3 < (eo p d) 2 /d 3 (perhaps tens of eV per molecule) where the Coulomb model's estimate 
(66) exhibits instead a nuclear Coulomb energy of tens or a few hundred MeV per nucleus (of which there are a few 
per molecule). Thus the true |V^| is a fraction < 10~ 7 (i?/r) 3 of the Coulomb model prediction. Hence, if it were 
possible to accurately measure T's acceleration at r ~ 10R from S, the ferroelectric would cause WEP violations 
at the fractional level < 10~ 16 (l/£p) 2 or smaller. This is within the projected sensitivity of the STEP experiments 
(which, however, will not carry aloft sources of tp). But because the electric field near a ferroelectric is very strong, 
on the order 10 8 (R/r) 3 V/cm for PbTi0 3 , it would probably prove necessary to work at considerably smaller R/r to 
avoid electric perturbations of T. The WEP violation would then most surely become unobservable at the projected 
STEP sensitivity. 

To sum up this subsection, there are no clear cases where extant or planned tests of the WEP are expected to 
detect violations connected with a variability which originate in the electric structure of matter. But a STEP-likc 
experiment based on a massive ferroelectric source of tp field in Earth's orbit would be a step closer for detection of 
WEP violation arising from a variability. 



C. Why does the Coulomb model fail so badly ? 

Why is the Coulomb model recapitulated in Sec. VA, with its seemingly unassailable logic, so inaccurate ? The 
answer to this puzzle hinges on a peculiar cancellation which parallels that which facilitated the passage from the 
complicated Eq. (26) to the simple Eq. (27) in the context of one charged particle. I now explain with the help of 
macroscopic averaging. 

In S each microscopic particle's mass is subject to Eq. (43). Therefore, macroscopic averaging of Eq. (38) for ip 
entails the replacement 



drfiiC 



uni-c 2 If 1 

X)— q -> -y / d 3 a;XK~ 1 eoitan(K$(z 4 ))<5 3 (x-z i ) w --^e ,*(z,) (68) 

where the approximation uses the fact (Sec. IIIC) that at the microscopic level n\<&\ <C 1, even near a charge, to 
discard terms of C(k 3 $ 3 ). Note that this expression is minus twice the formal microscopic Coulomb energy density 
in S; the sign is opposite that expected from the logic of the Coulomb model. 

Next one needs the macroscopic average of e~ 2 ^Y? in S. In the above mentioned approximation Eq. (46), E = 
-e 2 ^V$, gives with Eq. (45), e 2 ^ w 1 + k 2 $ 2 + • • • , that 



(69) 



^-e- 2 ^E 2 ^-±- { A[(V$) 2 l«¥E 2 ]4y%$W-A / QE-ds-n 2 [ d 3 x$ 2 E 2 
47T 47TV J v V t=v ^ dv ^ v 

The second form comes from integrating the (V$) 2 with the help of Gauss's theorem and employing Eq. (40). 

Note that the aforementioned microscopic Coulomb energy density term cancels out from the combined contributions 
(68) and (69). One term that survives is the surface integral over <i>E. To estimate it recall that if dV is pushed 
outward somewhat, $ on it becomes largely immune to fluctuations from individual charges composing S. Roughly, 
then, by Gauss's electric law the surface term in Eq. (69) is, V" 1 ($) J2iev e oi where ($) is an average of $ over the 
surface. But this expression is on the order of the Coulomb energy per unit volume associated with the net free charge 
contained in V. And we have seen by example in Sec. VB 2 that this last energy is much smaller than its microscopic 
counterpart V -1 J2iev e 0i$(zj)- 

What about the k 2 $ 2 E 2 integral in Eq. (69)? Let L be the smallest microscopic length scale on which $ varies, 
e.g. 10~ 13 cm if we think of nucleons as the smallest constituents, and perhaps 10~ 17 cm if quarks and electrons 
are considered in their stead (already in orthodox quantum electrodynamics, their electric fields are not scale-free 
Coulombic). Then one expects the maximum |$| to be no larger than a unit charge divided by L. Hence on the the 
average k 2 $ 2 E 2 /47t is bounded by a(l/L) 2 times the average microscopic Coulomb energy density E 2 /47r. This source 
of -0 is a fraction < 10~ 34 (l/£p) 2 of that assumed in the Coulomb model. It is thus amply clear how the cancellation 
of Coulomb energy source terms in the equation for -0 causes the Coulomb model to grievously overestimate |V-0|. 
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VI. MAGNETIC STRUCTURE AND WEP 



In ordinary matter, the next largest energy to Coulomb energy is energy of the magnetic dipoles associated with 
spin and orbital angular momentum of nuclei and electrons. Magnetic energy has so far been ignored here. Eq. (38) 
is missing the magnetic part of the term e~ 2 ^ f^f 1 *" in the original scalar equation (16). And in Sec. VC magnetic 
contributions have been left out of the mass term in Eq. (68). Now that it is clear that Coulomb contributions are 
well nigh irrelevant as sources of r/S, it is mandatory to take into account the purely magnetic contributions. 



A. Magnetic monopoles: a shortcut to dipoles 

It is a hard task to directly find the analog of the multicharge exact solution (45)-(46) for a collection of magnetic 
dipoles. But as the example of Sec. VC makes clear, this last solution is essential to derive the correct macroscopic 
form of the mass term in the source of Eq. (38) for ip. To overcome the problem I first find the analog of Eqs. (45)-(46) 
for a collection of magnetic monopoles of strengths goi, i — 1,2, •• • 2N , and then let pairs of equal and oppositely 
charged monopoles merge to form magnetic dipoles. The monopoles serve as a calculational crutch here and in the 
next subsection; they disappear from the final results. 

With the notation B = {/ 23 , / 31 , f 12 } I first replace Eqs. (37)-(38) by 

V x (er^B) = (70) 
V-B = 47T^ 5o ^ 3 (x-z 4 ) (71) 



V 2 -0 = 4ttk 2 



ydm^ 3 ( } _ J_ e -2^ B 2 
4^ dip V ' 4tt 



(72) 



Eq. (70) is the space-space component of the Maxwell-type Eq. (14) in a local Lorentz frame and with current and 
time derivatives set to zero. Eq. (71) docs not come from the action (7), but rather generalizes the Gauss law 
V • B = (which follows - also in variable a theory — from the representation of B as a curl) to the case that there 
are magnetic monopoles present. This is the argument for not including a factor e~ 2 ^ in (71). Finally, Eq. (72) is 
the scalar equation (16) in a local Lorentz frame with electric field and time derivatives dropped. 

It is not possible here, in light of Eq. (71), to use a vector potential for B; however, it can be written in terms of a 
scalar potential (its force on a monopole should be conservative): 

B = -V*. (73) 

It follows from Eq. (71) and the reasonable boundary condition W(x) — > as |x| — > oo that 

i 

Eq. (70) now gives e^ 2 ^"Vtp x = which uniquely implies that tp — ^>(\I/(x)). Finally, taking the Laplacian of 
this last relation gives 

v 2 ip = ip'vH + ip"(vy) 2 

= -47n//]T. 9ch <5 3 (x-z 4 )+^"B 2 (75) 

i 

Comparison of the second line of Eqs. (75) and Eq. (72) gives 

^P" = - K 2 e- 2 ^ (76) 
k 2 (dm,ic 2 1 dtp) = —goi?P' at x = z,. (77) 

The first of these integrates to 

V/ 2 = n 2 (e- 2 ^+m) (78) 
with w a constant. All the solutions of this equations may be obtained as in Appendix A: 



±( K * + x); w = o, 

±w- 1 / 2 smh.(w 1 / 2 K^ + x); w>0, (79) 

|W|- 1 / 2 C0S(|W| 1 /2 K $ + ^); ^< 
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with x a second integration constant. 

Because we made no guesses or ansatze along the way to Eqs. (79), these last together with Eq. (73) and (74) 
must exhaust the solutions of Eqs. (70)-(72) with boundary conditions ^(x) — > as |x| — ► 0. In this connection it is 
appropriate to remark that the substitution e 2 ^B E followed by ip *— > —ip and goi i— » e i, transform the magnetic 
solution (73)-(74) together with each of the five variants in Eq. (79) into one of the five variants of the electric 
solution (41)-(42) of Sec. IV A. A similar transformation maps each electric solution back onto a magnetic solution. 
The one-to-one correspondence between magnetic and electric solutions reflects duality symmetry of the a variability 
framework. Since all time-independent magnetic solutions have been found, the principle of duality assures us that 
all time-independent electric solutions have been found as well. 

Also by duality, the magnetic solution with w < 0, which is dual to the physical electric solution (w < 0), must 
be regarded as the physical choice. In fact, if invoking by now well worn arguments, we set w = — 1 and \ = 0, 
that solution is consistent with the reasonable boundary condition ip — > as |x| — > oo, whereas the solutions with 
vo > cannot be made so. Thus the unique physical solution for a collection of magnetic monopoles at rest is given 
by Eqs. (73) together with 

= cos(k#) (80) 
B = -V*. (81) 



B. Magnetic dipoles as source of ip 

Returning to our main subject, magnetic dipole energy, I imagine that the monopole labeled by i + N has magnetic 
charge —goi, and denote the vector from monopole i + N to monopole i by d. Then if as the pair are allowed to 
approach each other adiabatically, \g oi \ is made to grow in such a way that goid — > fi ^ in the limit, a point magnetic 
dipole is formed. It possesses magnetic moment fi and may be labeled by i. I still refer to its position as Zj. 

The total magnetic potential is 

go* -goi 

|x-z;-d/2| |x-z 4 + d/2|_ 

U |x-z,| 3 

Thus is explicitly known from Eq. (81). Performing the gradient in Eq. (81) gives for the field of the dipole, 

N 

b = E 

i=l 

This coincides (apart from the singularities at x = z$) with the Maxwellian expression for the field of an array 
of magnetic dipoles derived from the vector potential [28]. It seems reasonable that were magnetic dipoles to be 
represented by tiny current loops, the field would still be (83). Eqs. (80)-(82) thus represent the multimagnctic 
dipoles solution within the framework. 

What is the source term in the equation for ip for matter made exclusively of magnetic dipoles ? From Eq. (80) it 
follows that 

V 2 V> = k tan( K *)V • B - k 2 scc 2 (k*)B 2 (84) 

It is easy to verify, for example by applying Gauss' theorem to Eq. (83), that V • B vanishes everywhere, and has no 
5-function singularity at the position of a dipole. Since sec 2 (K\&) = e~ 2 ^, comparison of the last result with Eq. (72) 
shows that the mass-dependent source term is absent for magnetic dipoles. It is easy to get a wrong answer in this 
respect, for example, by trying to combine the mass-dependent terms of the positive and negative poles before taking 
the limit |d| — > 0. Since the mass-dependent source term would arise naturally from the action, the conclusion must 
be that the mass of a pointlike magnetic dipole, unlike that of a pointlike charge, does not depend on ip. 

An immediate consequence is that a magnetic dipole moment in an exterior electromagnetic field is not subject 
to an anomalous force of the sort appearing for an electric charge, cf. Eq. (26). The only force would be the usual 
(pi ■ V)B obtainable in our context by combining the forces on two magnetic poles. All this means that no violations 
of the WEP are expected in the motion of matter with pure magnetic dipole structure; it is unnecessary in this 
connection to carry out the analog of the calculation in Sees. IV B and IV C. 



lim > 
|dHof-f 



3/it • (x - Zj) (x - Zj) 



(83) 
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But does a source S with purely magnetic dipole structure cause WEP violations in the motion of a test body T 
made up of electric charges ? In analogy with the argument in Sec. VB it follows from Eq. (80) that 



1 



^,~__ K 2^2 = __ K 2^2__ K 2^_^2 ( 85 ) 

Contrary to what happens in the presence of a source with electric structure, cf. Sec. VB, here the physical a has 
its value slightly depressed in relation to the asymptotic value. As in Sec. VB, outside a macroscopic source S, the 
fluctuation term should be relatively small compared to 'J 2 , except in rather artificial situations where ^> vanishes (S 
has no net magnetic moment). The latter are not important in our context, so henceforth I shall drop the fluctuation 
and write just W for 

Supposing S to be a sphere with uniform magnetization M, ^ in Eq. (82) is well approximated at distance r >• R 
from the sphere's center by (47ri? 3 /3)M • rr~ 3 . Then 

|V^| ^4 K 2 (R/r) 3 [(4irR 3 /3)M 2 ]r- 2 (86) 

If Earth is the source <S, its crudely dipolar magnetic field (~ 0.25 G at the magnetic poles) allows its representation 
by a uniformly magnetized sphere with |M| ~ 3 x 10~ 2 cgs [28]. The energy represented by the square brackets in 
Eq. (86), <~ 9.7 x 10 23 erg, is a fraction ~ 10~ 22 of the corresponding factor in Eq. (66) for the Coulomb model (I take 
( = 10~ 3 for Earth). Hence Earth, by virtue of its magnetic structure, causes WEP violations at the fractional level 
10~ 2S (l/£) 2 , well beyond anything measurable in the foreseeable future. 

Another example of S is furnished by a ferromagnetic sphere. Iron, one of the ferromagnets with the highest 
saturation magnetization, can reach |M| w 1.7 x 10 3 cgs [30]. The energy in the square brackets in Eq. (86) is then 
s=s 1.2 x 10 7 erg/cm 3 R 3 . Being iron (p — 7.87 g/cm 3 and mass number 55.84) the sphere contains <~ 3.55 x 10 23 /cm 3 R 3 
Fe nuclei; each contributes <~ 125 MeV of Coulomb energy [29] for a total of ~ 7.11 x 10 19 erg/cm 3 R 3 . Thus the 
true IV^I, Eq. (86), amounts to a fraction 1.7 x 10~ 13 (i?/r) 3 of the Coulomb model prediction. It follows that a 
ferromagnetic source of \j) in orbit would cause WEP violations at a fractional level well below 10 -19 . This is beyond 
the projected STEP sensitivity (and, of course, the STEP experiment will not carry the equivalent of S). 

To be sure the above procedure is a bit cavalier. The exact solution (80)-(82) which forms the ultimate basis of the 
calculation contains no electric charges, while Earth and iron have plenty of them. But this defect is not serious. Let 
ip = ■0o( ( f ) ) + V'mW with ip m the function V'W m Eq. (80) and tp e the ip(&) defined by Eq. (45) for the collection of 
electric charges. Then the electric field E in Eq. (46) solves the Gauss equation (37), to within a relative correction 
of 0(V> m ), and is a gradient. Likewise, the magnetic field B [Eq. (83)] deriving from ^ in Eq. (82) satisfies, to within 
a fractional correction of (D(ip e ), the Ampere-like Eq. (70), and is also divergence free (Gauss's magnetic equation). 

Further, adding Eq. (38) for ip e to Eq. (72) (without the mass term) for ip m gives, to within fractional corrections 
of 0(tp c ) and 0(tp m ), the correct equation for the full tp relevant for a collection of charges and magnetic dipoles, cf. 
Eq. (16): 



V 2 ^ = 47r« 2 ^<5 3 (x - Zi ) + ^e- 2 ^(E 2 - B 2 )] (87) 



dip 47T 

The sum here extends only over electric charges. Plainly the electric and magnetic contributions to ip are additive 
(because \ip c \ <C 1 and |^ m | <C 1). Thus the correct V> to use for real matter (charges plus magnetic dipoles) is, not 
that in Eq. (85), but rather 

$ „ ($2 _ ^ (88) 

In view of this composition rule, the results of Sec. VI B that WEP violations originating from ip m will be unob- 
servable for the foreseeable future, together with those of Sec. VB that effects coming from ip e have the same status, 
show there are no clear cases where extant or planned tests of the WEP are expected to detect violations connected 
with a variability. Still outstanding is the question whether motion of charges or dipoles in matter — totally ignored 
in this paper could provide a loophole from this conclusion. 



VII. RESOLVING A COSMOLOGICAL CONUNDRUM 



The cosmological evolution of a according to the framework is influenced by the nature of ^'s source and by a 
certain integration constant tc (a cosmological time scale) [9]. If the source is described with the Coulomb model 
applied to ordinary matter of cosmological abundances, and the appealing choice tc — is made, then a is predicted 
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to decrease throughout the matter dominated era [9, 11, 14]. The claimed cosmological rise of a [1] thus led Olive 
and Pospclov [11] and Barrow, Sandvik and Magueijo [16] to conclude that the scalar field is coupled mostly to dark 
matter with a coupling opposite in sign to the naive one. To quote Ref. [16], "... the dark matter constituents have 
to have a high magnetostatic energy content (one possible contender would be superconducting cosmic strings which 
have C ~ —1)"; as clear from Eq. (87), magnetic fields are sources of tp with effectively negative (. 

Arguably the mentioned resolution of the cosmological conundrum amounts to trading one cosmic mystery (growing 
a) for another (dark matter with especially unusual electromagnetic properties). A much more natural solution is 
offered by the realization (see Sec. VC) that the Coulomb energy cancels from the source term of ^>'s equation. This 
should also be true for cosmological baryonic matter because one can think of each small region of the cosmological 
medium as being at rest in a local Lorentz frame. Consequently, the source of ijj's cosmological evolution is principally 
the B 2 term contributed by matter's magnetic dipoles to the cosmological version of Eq. (87). [By Sec. VI B the 
mass term in Eq. (16) is absent for magnetic dipoles, thermal radiation does not contribute at all because E 2 — B 2 
vanishes for electromagnetic radiation [9], and cosmological magnetic fields are too weak to make a difference.] Thus 
the switch in sign of the source term required by the claimed cosmological increase of a comes about automatically 
by considering the magnetostatic energy of baryonic cosmological matter. As I show now, the rate at which a grows 
also comes out of the right order without special assumptions. 

For cosmological a variability the interesting field is V>. As already discussed, it will obey an equation like (16) but 
with the source averaged over a macroscopic region. This average magnetic energy density is easiest evaluated by first 
computing the magnetic energy associated with a given mass of cosmological matter. 

Now in baryonic matter the magnetic energy is principally tied to the protons' and electrons' magnetic moments 
( 4 He nuclei have no magnetic moment). A lower bound on (and a passably good estimate for) the magnetic energy 
E p of the proton should result from regarding it as a perfect magnetic dipole and ignoring the energy interior to 
the proton radius R p . Integrating B 2 /87r (cf. the electric energy in Appendix B) with B in Eq. (83) over the space 
outside R p gives E p > /i 2 /3i? p 3 . Now for the proton /i = 2.98 x eh/2m p c [29], whereas R p w 5.7(h/m p c). Thus 
E p > 2.54 x 10~ 5 m p c 2 . For the electron /i = eh/2m e c. In view of the Zitterbewegung phenomenon is seems reasonable 
to integrate the magnetic energy of the dipole field only down to the Compton length h/m e c. Doing this as in the 
proton calculation gives E e > 6.1 x 10~ 4 m e c 2 . Other magnetic energies, e.g. those connected with spin-orbit and 
hyperfinc splittings, are small on the Bohr scale a 2 m e c 2 , and thus negligible here. 

In the context of a matter-dominated Robertson- Walker expanding model (with no cosmological constant and 
neglecting ip's effect on the expansion as appropriate at late cosmological time), one can use Eq. (31) of Ref. [9] for 
the rate of change of a: 

a/a = -(3/4tt)(//^p) 2 C n b H 2 t (89) 

Here H is the Hubble "constant" at cosmological time t, fit, the baryon density parameter and \(\ the fraction of the 
cosmic baryon mass density p b which acts as source of V> (here £ < 0). Now 23% of p b consists of 4 He. From the 
above results the magnetic energy density from helium's two electrons, and the protons and electrons from hydrogen, 
when appropriately weighted by their cosmic abundances, is S m ^ 1.98 x 10 ^ p b c 2 7 meaning that £ ^ — 1.98 x 10 ^. 
With Ht = 2/3 (matter domination) and today's favorite value VL b = 0.03, Eq. (89) predicts 

a/a > 6 x 10- 8 (//^p) 2 t _1 . (90) 

As advertised, a is predicted to increase with time. In order for its overall fractional change over the last 10 10 
y to match the observed 0.7 x 10 -5 , it is necessary that / < 10£p. This is a reasonable value consistent with the 
framework's assumptions and with all tests of the WEP (see Sees. VB and VI B). Hence nothing but baryonic matter 
is necessary to explain the claimed cosmological sense and rate of a variation. Were this last to be ruled out by future 
observations, one would have to conclude that I is smaller than mentioned. Within the spirit of the framework an I 
much below £p is unacceptable, and so one can risk a flat prediction that in recent epochs ta/a > 10~ 8 . Observational 
exclusion of this bound would be tantamount to certifying that standard Maxwellian electrodynamics (the case I = 0) 
is the exact classical description of electromagnetism. 

VIII. CONCLUSIONS AND CAVEATS 

The issue of a variability has long been connected with the possibility of weak equivalence principle violations 
coming from the classical Coulomb energy contribution to particle masses. The central conclusion of this paper is 
that in the general framework for a variability described in Ref. [9] , a compensating mechanism exists which prevents 
such violations from being measurable, at least in the foreseeable future. The treatment has been classical or tree 
level. Quantum considerations raise a potential problem for this kind of theory. 
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Banks, Dine and Douglas [31] have argued that in a theory where a scalar field couples to the electromagnetic 
scalar f^f^, as in the action (7), the residual vacuum energy of matter fields is very much larger than the observed 
cosmological constant unless the parameters are finely tuned. Recall that even with no scalar coupling the vacuum 
energy is known to be formally very large. But it is widely believed that some as yet ill understood mechanism nearly 
cancels it. Ref. [31] points out that because of the variation of the scalar field required for the suggested explanation 
of cosmological a variation to work, this cancellation can be successful only over a very short interval of cosmological 
time, with the residual vacuum energy becoming intolerably large earlier and thereafter. 

In light of this argument, can a theory like the framework used here be taken seriously ? One has to recall that the 
aforesaid criticism presupposes the existence of a cancellation of the vacuum energy by a mechanism whose nature is 
not generally agreed upon today. One can well imagine that when that mystery is lifted, a resolution for the problem 
of large variation of the residual vacuum energy with the scalar field might also become apparent. At any rate, it is 
premature to draw any final conclusions at this stage since the problem remarked upon in Ref. 31 is an integral part 
of the unsolved mystery of the small cosmological constant. 
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APPENDIX A: SOLUTIONS FOR V 

For w = take a square root of Eq. (24) to get 

e"V = ±« (Al) 

which may immediately be integrated to 

e^^TinV + x)- 1 (A2) 

with x a constant. 

For w ^ set tp = Iny + ^ In \zu\. Then after taking a square root, Eq. (24) transforms into 



dy 



V\Jy 2 + 1 • sgn(tn) 



±\w\ 1/2 ndV (A3) 



In the case w > the l.h.s. of this last equation is the differential of arccsch(y) [the upper (lower) signs corresponding 
to negative (positive) y\. Thus 

j/ = ±csch(tu 1/2 Ky + x)- (A4) 
In the case w < the l.h.s. of Eq. (A3) is the differential of arcsec(y) for either sign of y. Thus since sec(±y) — sec(y), 

r ±( K $ + x) _1 ; ^ = 0, 

= <; ± % /tncsch( % /n7K<I> + x); vo > 0, (A5) 

{ scc(vVI k$ + x); ™ < o. 



APPENDIX B: APPLICABILITY OF THE REISSNER-NORDSTROM METRIC 

Here I focus on an elementary charge. To assess the corrections to the metric coming from the scalar field I write 
down the energy-momentum tensor contributions from ip and f^ v in curved spacetime in accordance with Eqs. (6)-(7): 

%A = ^2 (t^- 1 -^,"^) (Bl) 
= -^ r [F a U--9^f al3 U) (B2) 
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In either tensor energy densities and stresses have the same magnitudes for a spherical static solution for which 
necessarily ip iT and f tr are the only surviving components of ip tli and f^ v , respectively). In the weak gravity region 
(r >• R) one may use Eq. (20), the fact that Vi/> = ^'V$ and Eq. (33) (with x = 0) to get for the said solution 

= «" V 2 e 2 ^ = sec 2 tan 2 (B3) 

T^^v is thus negligible compared to T^rfv all the way in to a radius where |$| is no longer small compared to 
7r 3 / 2 (?ic) 1 '' 2 Z _1 . According to Sec. IIIC the said radius represents the smallest possible extension of the charge. Hence, 
unless the charge is as compact as permitted, one can consistently neglect T^^ v in the above equations over the 
entire charge's exterior, and can do so over most of the exterior in the compact charge case. 

A look at Eq. (21) shows that ftr here is e 2 ^ times its Maxwellian counterpart cV<i>. Hence, T/f^v here differs 
from the pure Maxwellian energy momentum tensor also by a factor e 2 ^ = sec 2 which is very close to unity all 
the way in to the radius mentioned above. Thus the Maxwellian energy momentum tensor replaces the full T M „ 
throughout the whole exterior of a not maximally compact charge, and throughout most of it for a compact charge. 
As a consequence the Reissner-Nordstrom metric (35) is rather accurate outside the charge. 



APPENDIX C: THE MASS FUNCTION 



What can be said about spatial variation of masses of electrically charged particles through their dependence on ip 
? For simplicity of notation I focus on that of charge 1 in the cluster discussed in Sec. IV B. I first rewrite Eq. (43) 
with the help of Eq. (45) as 

(9mic 2 /a$) = -e i tan 2 («$). (CI) 

Evidently <& here means the value of the potential (39) evaluated at the position of charge 1, and of course this 
potential is large there (not infinite, though, since as mentioned in Sec. IV B, charges cannot be exact points in this 
framework). The equation refers to the change of m\ resulting from the local change of $ which might be due, for 
example, to the other charges being moved around (adiabatically) . 

Proceeding purely formally I use the fact that dt&nx/dx = 1 + tan 2 x to integrate Eq. (CI): 

mic 2 = m ic 2 + e i [$ - KT 1 tan(/c$)] (C2) 

Here moi is a constant of integration. It is clear from this that mi is invariant under conjugation of all charges, as 
would be expected by C invariance of the electromagnetic interaction; although the framework's assumptions require 
only T invariance, the theory's action (6)-(8) is actually C invariant as well. According to Eq. (33), n<& = arcsece^; 
since tan 2 (K$) = sec 2 (ft$) — 1 one has 

mic 2 = m ic 2 + K _1 e i ^arcsece^ — \J e 2 ^ — lj . (C3) 

The appropriateness of choosing the shown sign of the square root is verified by differentiating mi (tjj) and comparing 
the result with Eqs. (43) and (33). The opposite sign is excluded on this ground. 

As anticipated, the space-dependent contribution to mi is not proportional to e i 2 e 2 ^ as would be assumed in the 
Coulomb model. 

Unfortunately, Eqs. (C2) and (C3) are not immediately useful; they contain a contribution from the self potential 
of particle 1 or from its self ip field. If there were no other charges in the universe, one would expect mi to subsume 
these contributions. Absorbing the self terms of Eq. (C2) into moi gives 

mic 2 = m ic 2 + eoi{$i>i ~ [tan(ft$) - tan(/c$j=i)] } (C4) 

The identity (tana; — tany) cos £ cosy = sin(x — y) allows us to rewrite this as 

mic 2 = m ic 2 + eoi[$i>i - n~ x A sin(K<I> l> i)] (C5) 

where A = sec(«<I> i= i) sec(K<I>) is evaluated at particle l's center. 

Quite generically k|$,>i| <SC 1. This merely requires that the overall size R of the cluster of charges i — 2, 3, • • • , 
and its overall charge sans eoi, Q, satisfy 



R^7T- 3 / 2 (Q 2 /hc) 1 ^(l/£ P )e P 



(C6) 
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which is easily met by nucleons, nuclei, ions and macroscopic charged objects. One expects charge 1 to be relatively 
more compact so that |$i=i| S> |$»>i|- In this case Taylor expanding A to first order in K$i>i as well as sin(«;$ i>1 ) 
to second order converts Eq. (C5) into 



mic 2 = m 01 c 2 - e 01 scc 2 (k$ 4=1 ) x [sin 2 ^^)^! + /ttan(/«$ i=1 )$ i>1 2 + O^i^ 3 )} (C7) 

With the possible exclusion of the "point" leptons and quarks, even elementary particles are much more extended than 
the bound (34); thus fc|<I>j = i| <C 1. Hence the variable part of mi is generally very small compared to the Coulomb 
potential energy of that charge in the cluster, eoi^^i, and of opposite sign. 
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